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摘要 : 锥 束 CT 的 图 像 阴影 问题 严重 影响 了 其 CT 值 精度 、 低 对 比 度 目 标的 检测 能 力 及 剂量 计算 的 准确 性 
限制 了 锥 束 CT 在 临床 上 广泛 应 用 。 为 了 消除 图 像 阴影 ， 文 章 首 先 对 原始 CT 图 像 进行 多 阔 值 的 图 像 分 
害 ， 把 骨骼 和 软组织 分 离 并 赋予 标准 的 CT 值得 到 模版 图 像 ， 然后， 利用 Lo 范 数 平滑 算法 对 原始 CT 图 
像 进行 边缘 保护 平滑 ， 得 到 去 除 图 像 纹理 信息 的 平滑 图 像 ， 之 后 两 幅 图 像 做 差 后 利用 低 通 滤波 即 可 得 到 图 
像 阴影 分 布 ， 最 后 把 估计 的 阴影 分 布 补偿 至 原始 图 像 ， 得 到 阴影 校正 图 像 。 实 验 发 现 校正 后 的 感 兴趣 区 域 
的 CT 误差 由 大 于 115 HU 下 降 至 小 于 13 HU， 整 体 图 像 非 均匀 度 由 大 于 996 下 降 至 小 于 1%。 实 验 结 
RRI, EF Lo 范 数 平滑 和 图 像 分 割 的 CT 图 像 阴影 校正 方法 可 以 有 效 地 校正 CT 图 像 阴影 ， 具 有 一 定 
的 临床 应 用 价值 。 
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Abstract: X-ray shading artifacts lead to CT number inaccuracy, image contrast loss and spatial non-uniformity, and therefore are 
considered as one of the fundamental limitations of CBCT. In order to solve this problem, a novel shading correction method was 
proposed. Specifically, we first use multi-threshold segmentation algorithm to segment the original CT image for constructing a 
template image where each structure is filled with the same CT number of a specific tissue type. Then, the Lo norm smoothing 
algorithm is used to smooth the CBCT image for constructing an image without texture. By subtracting the template from the image 


without texture, the residual images from various error sources are low-pass filtered to generate the estimated shading artifacts. 
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Finally, the estimated shading artifacts are added back to the original image for shading correction. Compared with the CT image 
without correction, the proposed method reduces the overall CT number error from over 113 HU to be less than 13 HU and decreases 
the non-uniformity from over 996 to be less than 196. The experimental results demonstrate that the proposed shading correction 
method using Lo norm smoothing and image segmentation can effectively correct the shading artifacts and its feasibility in clinical 
application is validated. 
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1 引言 
在 锥 形 束 CT (Cone-Beam CT, CBCT) 扫描 中 ， 由 于 散射 信号 由 和 射 束 硬化 效应 中 的 影响 ， 重 建 图 像 
中 出 现 低频 的 CT 图 像 阴影 ， 这 严重 影响 了 图 像 CT 值 的 准确 性 和 图 像 的 空间 均匀 度 S9。 在 无 阴影 校正 
的 CBCT 系统 中 ，CT 图 像 阴影 导致 重建 图 像 的 CT 值 误差 可 超过 350 HU 中， 给 影像 引导 治疗 的 定位 精 
度 和 影像 的 诊断 带 来 误差 ， 因 而 限制 了 CBCT 在 临床 中 的 广泛 应 用 。 目 前 CBCT 主要 还 是 用 于 初步 定位 
和 放疗 摆 位 ， 在 介入 和 放疗 中 的 进一步 应 用 受到 限制 ， 因 此 CT 图 像 阴影 校正 是 提升 CBCT 图 像 质量 首 
要 解决 的 重要 问题 之 一 。 

目前 已 知 的 CT 图 像 阴影 校正 的 方法 主要 分 为 两 大 类 : 预 处 理 和 后 处 理 方法 四。 其 中 ， 预 处 理 方法 校 
E CT 图 像 阴影 主要 是 通过 附加 硬件 装置 ， 以 阻止 散射 光子 到 达 探 测 器 ， 这 样 也 就 不 存在 散射 信号 。 预 处 
理 进行 阴影 校正 有 两 个 典型 方法 : 一 是 增加 物体 与 探测 器 之 间 的 空气 阶 四 ， 二 是 使 用 抗 散射 线 栅 昌 。 随 着 空 
CAUSE. 扩散 开 的 散射 光子 的 探测 率 会 降低 ,而 源 信号 则 不 受 影响 。 但 该 方法 受 CBCT 设备 本 身 物理 
间 的 限制 ， 空 间距 离 不 能 无 限 增 大 ， 而 且 增 大 的 同时 加 大 了 影像 的 几何 模糊 ， 同 时 还 需要 增加 X 线 剂 量 
来 弥补 距离 的 增加 ， 在 临床 实际 中 并 不 实用 。 抗 散射 线 栅 使 用 聚焦 于 射线 源 的 铅 栅 网 格 ， 能 阻挡 非 聚 焦 入 
射 角 的 散射 光 。 该 方法 也 存在 对 散射 光 的 衰减 效率 不 高 的 缺陷 。 目 前 商用 铅 栅 只 能 提供 约 3 倍 的 散 透 比 
(Scatter-to-Primary Ratio, SPR) 降低 率 ， 无 法 保证 高 散射 环境 下 的 CBCT 图 像 质量 。 此 外 ， 它 还 需要 
增加 病人 的 受 照 剂量 来 补偿 被 衰减 掉 的 源 射线 强度 ， 临 床 应 用 价值 不 高 。 
虽然 预 处 理 方法 能 直接 阻止 散射 光子 到 达 探 测 器 ， 但 其 局 限 性 比较 突出 ， 相 比 而 言 ， 后 处 理 方法 更 实 
用 。 后 处 理 指 的 是 按照 原来 的 方法 采集 X 射线 投影 后 ， 再 通过 后 处 理 估计 散射 分 布 ， 通 过 源 投影 减 去 估 
计 出 的 散射 分 布 ， 即 可 对 其 进行 阴影 校正 。 后 处 理 方法 包括 : 解析 建 模 法 钙 、 蒙 特 卡 罗 模 拟 法 040、 源 调 
制 法 123、 测量 法 、 基 于 先 验 数 据 校正 法 、 基 于 极 坐标 的 校正 法 和 自 适 应 友 代 阴影 校正 法 。 其 中 ， 解 析 建 模 
法 外 认为 散射 信号 是 源 信号 通过 散射 核 后 的 响应 ， 计 算 速 度 很 快 ， 但 相应 的 散射 估计 精度 有 限 且 对 复杂 物 
体 需要 繁琐 的 调整 参数 ， 蒙 特 卡 罗 模 拟 法 425 是 散射 估计 的 “ 金 标 准 ”， 但 该 方法 计算 量 极 大 ， 十 分 耗 时 ; 
源 调制 法 0 是 在 X 射线 源 和 物体 之 间 加 入 高 频 调 制 器 ， 根 据 散 射 和 源 信 号 不 同 的 响应 特性 ， 在 频 域 上 把 
它们 分 离开 ， 但 该 方法 对 调制 板 的 制造 精度 要 求 高 ， 其 临床 应 用 效果 受到 实际 物理 因素 的 限制 ， 测 量 法 是 
在 X 射线 源 和 受 照 物 之 间 插 入 源 射 线 阻挡 器 (通常 是 铅 条 )0311, 这 样 探 测 器 上 形成 仅 包含 散射 信号 的 阴影 
区 ， 但 该 方法 要 改变 系统 的 硬件 设置 ， 操 作 难 度 较 大 ， 基 于 先 验 数据 校正 法 虽 能 较 好 地 得 到 校正 图 像 08， 
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该 方法 需要 对 


昔 助 放射 治疗 中 额外 的 ? 


E 验 病 患 信息 ， 因 此 无 法 作为 大 体积 CT 成 像 系统 阴影 校正 
决 方案 ; 基于 极 坐标 的 校正 法 通过 CT 图 像 阴影 在 极 坐标 下 的 分 布 特征 09,， 估 计 出 CT 
图 像 进行 极 坐标 转化 和 捐 


i 值 操作 ， 耗 时 长 ， 自 适应 迭代 阴影 校正 法 无 需 先 验 


该 方法 需要 对 重建 图 像 进行 反复 正 投影 操作 ， 计 算 效 率 低 。 
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的 方法 ， 提 出 一 种 新 的 CT 图 像 阴影 校 


FE 方法。 该 方法 不 


依赖 其 他 影像 数据 ， 完 全 兼容 现代 放疗 加 速 器 的 图 像 引 导 系 统 ， 不 用 改变 其 他 硬件 和 扫描 协议 ， 无 需 进 行 
坐标 变换 和 正 投 影 ， 不 增加 扫描 时 间 和 扫描 剂量 。 通 过 Lo 范 数 平滑 分 离 图 像 纹理 并 保护 图 像 细节 信息 ， 
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24 CT 图 像 阴影 校正 方法 的 流程 


图 1 显示 的 是 CT 图 像 阴影 校 


图 像 ， 估 计 出 CT 图 像 阴影 


的 分 布 ， 最 终 达 到 CT 


图 像 阴影 校正 的 目的 。 


FE 方法 的 工作 流程 图 。CT 图 像 阴影 信号 属于 低频 信号 ，CT 


图 像 的 细 


节 信 息 同样 属于 低频 信号 。 为 了 保护 图 像 细 节 信 息 以 及 更 精确 的 估计 CT 图 像 阴影 的 分 布 ， 本 研究 首先 利 


用 Lo 范 数 平滑 ， 在 对 图 像 进 
类 型 的 组 织 的 CT 值 大 小 基本 一 致 。 根 和 
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像 。 参 考 


构 误差 。 


误差 去 除 ， 从 而 得 到 低频 的 CT 
像 的 生成 以 及 低 通 滤波 。 
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HX — BRA 
图 像 所 含 结构 与 无 纹理 的 平滑 图 像 做 差 的 残留 
因为 结构 误差 属于 高 频 信 号 ， 而 CT 图 像 阴影 属于 低频 信号 ， 因 
影 分 布 。 下 面 将 对 工作 流程 中 的 关键 步 又 进行 介绍 ， 包 括 : Lo 范 数 
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图 1 图 像 阴影 校正 方法 的 工作 流程 图 
Fig. 1 Workflow of the proposed method of CT shading correction 


图 像 存在 CT 


uUum 识 指 出 , 


图 像 ， 模 版 


BEEN ERG 


EP Ne 


人 体 同 种 


图 像 可 以 作 
图 像 阴影 和 
此 可 通过 低 通 滤波 把 少量 的 结构 


少量 的 结 


22 Lo 范 数 平滑 

在 图 像 阴影 校正 过 程 中 ， 为 了 保护 图 像 细节 ， 不 损失 图 像 空间 分 辩 率 ， 本 研究 把 CT 图 像 分 解 为 纹理 
和 结构 。 相 对 于 CT 图 像 阴影 信号 而 言 ， 图 像 的 纹理 属于 高 频 信 号 。 为 了 更 精确 地 提取 散射 信号 ， 本 研究 
首先 对 原始 重建 图 像 进行 边缘 保护 的 Lo 范 数 平滑 。 有 具有 边缘 保护 的 Lo 范 数 平滑 算法 过 程 如 下 : Wt I 为 输 
0,S,) 为 图 像 的 梯度 Cp. 为 图 像 的 像 


入 的 原始 CT 重建 图 像 ，S 为 输出 无 纹理 的 平滑 图 像 ，V5S, 二 (0、S 


p? 


RRI x 和 y 分 别 为 图 像 的 横 纵 两 个 方向 的 坐标 )。 对 于 一 个 像素 索引 p. 计算 其 在 x 和 oy 两 个 方 
向 的 梯度 差 之 和 作为 它 的 梯度 值 ， 如 公式 〈1) PUB: 


Cc(s)=#{p| |ess,| 上 |eySn| 关 0 «o 


其 中 ，C(S ) 为 


6、S,| 十 |6, Si,| 产 0 像素 索引 p 的 个 数 。 根 据 上 面 的 定义 ， 无 纹理 的 平滑 图 像 S 可 以 通过 


以 下 目标 函数 求 得 C0， 见 公式 (0): 


vin [Y(5,-1 +2.c()| (2) 
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Hp, Y(s—rIy 为 图 像 结构 相似 性 约束 项 。 因 Lo 范 数 的 不 可 解 性 ， 所 以 本 文 将 算法 分 为 两 个 子 问题 进 


=>% 
Hl 


分 别 求解 ， 采 用 半 二 次 分 裂 的 特殊 交 葵 优化 策略 得 到 近似 的 最 优 解 。 


2.3 模版 图 像 的 生成 

根据 解剖 常识 ， 同 种 类 型 的 人 体 组 织 在 没有 伪 影 的 CBCT 图 像 中 的 CT 值 应 该 基本 一 致 上 四。 由 此 ， 可 
以 利用 图 像 分 割 的 方法 ,把 人 体 不 同 的 组 织 分 割 出 来 。CT 图 像 中 的 结构 成 分 可 分 为 空气 、 骨头、 软组织 (如 
肌肉 和 脂肪 ) 区 域 等 ， 对 应 的 组 织 区 域 填 充 相 应 组 织 在 对 应 Xx 射线 球 管 电压 下 的 标准 CT 值 ， 生 成 一 幅 模 
版 图 像 ， 该 图 像 作 为 阴影 校正 过 程 的 参考 标准 。 本 研究 采用 多 阔 值 的 图 像 分 割 算 法 ， 分 离 病人 头 部 的 骨骼 和 
软组织 并 赋予 标准 的 CT 值 。 


2.4 低 通 滤波 

上 述 生 成 的 模板 图 像 与 无 纹理 的 平滑 图 像 之 差 包 含 图 像 阴影 和 少量 组 织 结构 误差 。 阴影 伪 影 主要 为 低频 
信号 ， 而 组 织 结构 则 主要 是 高 频 信 号 ， 可 利用 低 通 滤波 器 消除 组 织 结构 误差 。 普 通 的 低 通 滤波 可 以 消除 组 织 
结构 的 误差 ， 但 解剖 结构 的 边界 轮廓 也 会 被 滤 除 ， 导 致 图 像 对 比 度 的 严重 丢失 。 因 此 ， 本 研究 采用 
Savitzky-Golay 局 部 低 通 滤波 器 22 在 图 像 域 上 对 残 差 图 像 进行 滤波 ， 该 低 通 滤波 器 可 以 保持 轮廓 特征 ， 避 免 
了 校正 后 损失 图 像 的 对 比 度 分 辩 率 。 


2.5 实验 评估 

通过 Varian Trilogy 机 载 CBCT 系统 采集 实验 所 需 投 影 ， 在 均匀 模 体 上 进行 CT 图 像 阴影 校正 ， 再 
对 病人 头 部 重建 图 像 进行 CT 图 像 阴影 校正 ， 最 后 ， 把 扇 束 CT 和 计划 CT 扫描 图 像 作为 均匀 模 体 和 病 
人 头 部 CT 图 像 阴影 校正 效果 的 金 标准 与 校正 后 的 图 像 进行 对 比分 析 。 对 于 扇 束 CT， 本 研究 把 射线 源 前 
放 的 准 直 器 宽度 缩小 为 4mm， 采集 均匀 模 体 ， 得 到 一 组 近似 无 CT 图 像 阴影 的 扇 束 CT 图 像 ， 即 可 利用 
AR CT 图 像 与 校正 后 的 图 像 作对 比 。 此 外 ， 实 验 采集 的 Capthan 504 模 体 数据 用 于 算法 参数 设置 分 析 。 
实验 采集 的 三 组 数据 扫描 与 图 像 重建 参数 如 表 1 所 示 。 图 像 处 理 的 硬件 配置 如 下 : Corei7-6820HK E7200 


2.70 GHz CPU、32 G 内 存 的 微星 笔记 本 电脑 ， 编 程 环境 采用 Matlab R2016a. 


ng 


ANS 


表 1 瓦 里 安 Trilogy 机 载 CBCT 系统 扫描 与 图 像 重 建 参 数 


Table 1 Imaging and reconstruction parameters of the Varian Trilogy OBI system 


Z3 均匀 模 体 Catphan 504 病人 头 部 

扫描 模式 full-fan full-fan full-fan 

= X 射线 能 量 (kVp) 125 125 100 

= 管 电流 (mA) 80 80 80 
X 射线 焦 斑 尺寸 (Omm) 0.4 0.4 0.4 

- 脉冲 宽度 Cms) 25 25 13 

-— 源 至 探测 器 距离 (mm) 1 501.42 1 501.42 1501.42 
源 至 等 中 心 距离 (mm) 1 000 1 000 1 000 
探测 器 尺寸 400*300mm? 400*300mm? 400*300mm? 
探测 器 像素 1024*768 1024*768 1024*768 
旋转 角度 C ) 360 360 360 

投影 数量 663 662 367 

重建 体 素 大 小 (mm Xmmx mm) 0.5X0.5X 0.5 0.50.5 X 0.5 0.50.5 X 0.5 
重建 体积 大 小 〈 体 素 ) 512X 512X 512 512X 512X 512 512X 512X 512 
算法 FDK FDK FDK 


在 图 像 的 定量 分 析 中 , 采用 CT 值 误差 与 整体 图 像 非 均 匀 度 评价 校正 前 后 图 像 质量 。 整 体 图 像 非 均 匀 


H£ (Spatial Non-Uniformity，SNU) 的 计算 公式 075 如 公式 GB) 所 示 。 


sNU — HU su HU vin y 40095. (3) 
1000 


在 图 像 同一 组 织 区 域 选 择 5 个 半径 大 小 为 5 个 像素 的 感 兴趣 区 域 CRegion of Interes; ROD. t, 


4 个 ROI 选 在 重建 物体 的 边缘 ，1 个 ROI 选 在 重建 物体 的 中 心 。 式 (3) 中 HU sac 和 HUS 分 别 为 


5 个 ROI 中 平均 CT 值 的 最 大 值 和 最 小 值 。 
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3.1 均匀 模 体 


对 比 图 。 从 图 2 和 图 


T 


图 2 为 均匀 模 体 中 心 层 的 横断 图 像 ， 图 3 为 阴影 校正 前 后 在 图 2 Co) 中 白色 虚线 位 置 上 CT 值 的 


3 可 看 出 ， 图 像 阴影 明显 地 被 抑制 ， 抑 制 后 的 图 像 效 果 图 2 CO 接近 于 扇 束 CT 


limi 
um 


建 的 图 像 。 图 2 (a) 中 5 个 白色 虚线 圆圈 为 计算 整体 图 像 SNU 的 ROIL 校正 后 图 像 的 整体 SNU 由 996 


E 0.4%; 图 2 Cc) 


EERE ROI 的 CT 值 误 差 从 77 HU 降 至 2HU。 此 外 , 图 2 中 3 个 图 像 的 


ROI (白色 虚线 圆圈 内 ) 的 CT 平均 值 分 别 为 : (a)—139 HU; (b) 一 60 HU; (QO 一 62 HU. E 3 表明 校正 后 
CT 值 强度 基本 吻合 。 


的 CT 值 强 度 与 扇 束 


(a) 无 CT 图 像 阴影 校正 的 CBCT (b) 校正 后 的 CBCT (c 参考 扇 束 CT 
图 2 均匀 模 体重 建 的 横断 面 图 像 〈 显 示 窗 口 为 : [一 150 450] HU) 
Fig.2 Axial views of the reconstructed uniform acrylic cylindrical phantom (Display window: [一 150 450] HU) 


CT 人 (HU) 


i 一 一 参考 户 来 CT 图 像 
seven 未 校正 的 CT 图 像 
-200r 一 一 校正 后 的 CT 图 像 1 


i i — i] 
0 50 100 150 200 250 300 


3 CT 图 像 阴影 校正 前 后 通过 图 2(o) 白 色 虚 线 所 在 位 置 的 CT 值 对 比 图 


Fig. 3 


3.2 TR A SERI 


Comparison of 1D profiles passing through the dotted white line as indicated in figure 2(c) 


图 4 显示 的 是 病人 头 部 的 横断 图 像 ， 图 5 Ca) 和 (b) 分 别 显示 阴影 校正 前 后 在 图 4 Cc) WU (D 中 


白色 虚线 位 置 上 CT 值 的 对 比 图 。 从 图 4 和 图 5 可 看 出 ， 图 像 阴影 明显 地 被 抑制 ， 抑 制 后 的 图 像 效 果 图 


4 第 二 列 接近 于 配 准 后 的 计划 CT 重建 图 像 。 图 4 GO 中 5 个 白色 虚线 圆圈 为 计算 整体 图 像 〈 图 a~ 


HRM 200 HU KÆ 13 HU。 此 外 ， 整 体 图 像 
173 HU; (c) 160 HU。 图 5 显示 校正 后 的 CT 值 强度 接近 于 配 准 后 的 计划 CT 值 强度 。 


图 4(O) SNU 的 ROI， 校 正 后 图 像 的 整体 SNU 由 1496 降 至 1%, 图 4 (b) 白色 虚线 圈 ROI 的 CT 
的 3 个 图 像 的 CT 值 平均 值 分 别 为 : (a) —27 HU; (b) 


(a) 校正 前 《中心 层面 ) 图 像 ; O) 校正 后 〈 中 心 层 面 ) 图 像 ; (oO 配 准 后 的 参考 计划 CT LED) Až; (0 校正 前 〈 非 中 心 层 面 ) 图 像 ; 


Ce) 校正 后 〈 非 中 心 层面 ) 图 像 ; (D 配 准 后 的 参考 计划 CT 〈 非 中 心 层面 ) 图 像 
4 病人 头 部 横断 面 CT 重建 图 像 〈 显 示 窗 口 为 : [一 200 500] HU) 
Fig.4 Axial views of the reconstructed patient head (Display window: [一 200 500] HU) 
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(a) 通过 图 4(O 白 色 虚 线 所 在 位 置 的 CT 值 (b) 通过 图 4(9) 和 白色 虚线 所 在 位 置 的 CT 值 
5 ”对 比 图 4 两 条 虚线 所 在 位 置 的 一 维 轮廓 像素 值 
Fig.5 Comparison of 1D profiles passing through the two dotted white lines as indicated in figure 4 
4 讨论 


Lo 范 数 平滑 的 参数 选择 影响 图 像 纹 理 去 除 程度 与 组 织 结构 边缘 保护 的 效果 , 纹理 去 除 和 边缘 保护 


对 矛盾 。 充 分 的 图 像 纹 理 去 除 ， 有 助 于 校正 后 图 像 的 细节 保 
助 于 估计 图 像 阴影 的 准确 性 ， 提 高 校正 后 图 像 组 织 间 的 对 比 


因子 4 的 选择 应 在 两 者 间 取 得 平衡 。 不 同 4 选择 对 图 像 造成 的 影响 如 图 
(cO 分 别 为 3 个 不 同 4 下 的 Lo 范 数 平滑 图 像 ， 图 6 (D —6 CO 为 在 相应 4 下 的 校正 


n 


导致 组 织 结构 边缘 模糊 ， 估 计 出 的 阴影 信息 中 包含 结构 信息 


当 A HX 0.01 时 ， 在 纹理 去 除 和 边缘 保护 之 间 取 得 了 平衡 ， 


度 。 


留 ， 不 破坏 图 像 分 辩 率 。 良 好 的 边缘 保护 ， 有 


平滑 强度 


取得 过 小 时 ， 如 图 6 (a) 和 6 (d)， 纹 理 去 除 不 充分 导致 估计 出 的 阴影 信号 中 包含 图 像 
从 至 原始 图 像 得 到 的 校正 图 像 损 失 图 像 细 节 。 当 4 取得 过 大 时 ， 如 图 6 (oO 和 65. E 
， 因 此 校正 后 图 像 组 织 结构 之 i 


, 


度 。 因 此 ， 在 Lo 范 数 平滑 过 程 中 ， 
6 所 示 。 其 ! 


图 6 (a) 一 6 


图 像 。 当 A 


节 ， 阴 影 信号 补 


平滑 程度 过 强 
司 对 比 度 降 低 。 


在 保护 图 像 细 节 的 同时 ， 不 降低 组 织 间 对 比 


(a) 平滑 图 像 (CA —0.00); ©) 平滑 图 像 (A —0.00; CO) 平滑 图 像 CA =0.04); (d) 校正 后 图 像 (4 三 0.001); (e) dz 
正 后 图 像 CA-000D; (D 校正 后 图 像 CA-0.04) 


6 不 同 平滑 强度 因子 4 下 的 平滑 和 校正 图 像 〈 显 示 窗 口 为 : [一 200 500] HU) 


Fig.6 Effect of the strength of smoothing A in the smoothing images and corrected images (Display window: [一 200 500] 


HU) 


精确 模版 图 像 的 生成 依赖 于 图 像 分 割 的 精度 。 本 文采 用 基 


T BA Otsu 准则 的 阔 值 分 割 


方法 ， 阔 值 由 最 大 类 间 准 则 确定 。 在 分 割 步骤 中 ， 需 要 假设 图 像 中 存在 n 个 待 区 分 的 类 ，m 
的 取 值 影像 图 像 阴影 修正 效果 。 在 病人 头 部 的 阴影 校正 中 , n 取 3 能 准确 分 害 
和 空气 ， 获 得 精确 的 模版 图 像 。 当 n 过 小 时 ， 无 法 区 分 灰 度 值 相近 的 不 同 组 


正 后 图 像 的 组 织 间 对 比 度 下 降 。 


阴影 。 与 类 似 研 究 相 比 ， 如 Niu 等 8 引入 诊断 CT 数据 ， 


过 大 时 ， 则 导致 同一 


织 


1 骨骼 、 软 组 


pmi 


， 导 致 阴影 修 


组 织 内 阴影 校正 不 全 ， 残 留 阴影 
本 研究 提出 了 一 种 结合 Lo。 范 数 平滑 与 图 像 分 割 相 结合 的 CT 图 像 阴影 校正 方法 ， 有 效 地 消除 了 图 像 


把 诊断 CT 作为 参考 图 像 对 CBCT 进行 阴影 


(S 


校正 ， 模 体感 兴趣 区 域 的 CT 值 误差 降 至 17 HU。 但 该 方法 对 诊断 CT 数据 的 依赖 很 强 ， 专 用 于 放疗 中 
的 影像 引导 , 并 且 由 于 诊断 CT 与 锥 束 CT 是 在 治疗 的 不 同 阶段 扫描 得 到 , 所 以 需要 对 二 者 进行 图 像 配 准 ， 
对 配 准 精度 要 求 高 。 而 在 本 研究 中 ， 模 体感 兴趣 区 域 的 CT 值 误 差 降 至 4HU， 主 要 原因 在 于 无 需 配 准 ， 
避免 了 配 准 误差 所 带 来 的 阴影 估计 误差 。 在 Fan 等 9 的 研究 当中 ， 首 先 把 原始 CBCT 图 像 由 笛 卡 尔 坐 标 
系 转 为 极 坐 标 系 ， 之 后 在 极 坐标 图 像 上 进行 滤波 ， 估 计 出 图 像 阴 影 ， 最 后 把 估计 的 图 像 明 影 转化 为 笛 卡 尔 
坐标 补偿 至 原始 CBCT 图 像 中 。 在 校正 后 的 病人 头 部 图 像 中 ， 同 一 组 织 的 图 像 SNU FEE 3%, 但 是 该 
方法 对 于 组 织 结构 复杂 以 及 带 环 形 的 物体 效果 不 好 .图 8 为 上 述 算法 和 本 文 提出 算法 的 Catphan 504 模 体 
处 理 结果 对 比 图 。 在 极 坐标 系 下 处 理 上 述 图 像 会 引入 新 的 伪 影 〈 如 图 8 (co 红色 箭头 所 指 )， 导 致 图 像 对 
比 度 的 损失 。 本 算法 处 理 的 效果 如 图 8 (b) 所 示 ， 保 持 了 图 像 的 对 比 度 ， 主 要 原因 在 于 引入 图 像 分 割 产生 
了 精确 的 模版 图 像 ， 模 版 图 像 的 同一 组 织 的 图 像 SNU 为 0， 再 利用 Lo 范 数 平滑 保留 图 像 细 节 ， 较 好 地 
保护 了 图 像 对 比 度 。 


(a) 无 CT 图 像 阴影 校正 的 CBCT (b) 本 文 提出 算法 的 校正 CBCT (c) 基于 极 坐标 系 算法 的 校正 CBCT 红色 得 
头 所 指 为 极 坐标 系 下 处 理 上 述 图 像 引 入 的 伪 影 ) 
图 8 不 同 阴影 校正 算法 的 效果 对 比 图 〈 显 示 窗口 为 : [一 30 450] HU) 
Fig.8 The results of the shading corrected CBCT using different algorithm (Display window: [一 30 450] HU) 


本 方法 利用 边缘 保护 的 Lo 范 数 平 滑 算法 ， 把 图 像 分 解 为 结构 图 像 和 纹理 图 像 ， 把 纹理 信息 消除 后 再 
做 后 续 处 理 ， 不 损失 图 像 分 辨 率 。 精 确 的 图 像 分 割 算 法 成 功 地 将 人 体 组 织 进行 分 割 ， 生 成 精确 的 参考 图 像 。 
因此 ， 可 作为 CBCT 图 像 阴影 校正 的 通用 方法 。 基 于 Lo 范 数 平滑 和 图 像 分 割 的 CT 图 像 阴影 校正 技术 
除了 图 像 阴影 校正 效果 明显 外 ， 还 具有 以 下 优点 : CL 通过 Lo 范 数 平滑 最 大 限度 地 保护 了 图 像 的 细节 信 
息 ， 图 像 空 间 分 辩 率 损失 小 ; (20 本 研究 无 需 先 验 的 CT 信息 ， 无 需 前 投影 操作 ， 无 需 极 坐标 转换 ， 因 此 
计算 速度 较 快 ; (3) 完全 兼容 现代 放疗 加 速 器 的 图 像 引 导 CBCT 系统 ， 不 用 改变 其 他 硬件 和 扫描 协议 。 在 
未 来 ， 本 方法 还 需 在 以 下 两 方面 进行 改进 。 第 一 ， 精 确 模版 图 像 的 生成 依赖 图 像 分 割 的 精度 ， 分 割 精度 更 
高 的 算法 有 助 于 图 像 阴影 的 校正 。 多 阔 值 图 像 分 割 算法 成 功 地 分 割 出 了 头 部 骨骼 、 软 组 织 和 空气 ， 但 难以 
分 割 出 脂肪 ， 导 致 校正 后 图 像 的 脂肪 和 其 他 软组织 的 图 像 对 比 度 下 降 。 因 此 ， 在 之 后 的 工作 将 进一步 提升 
图 像 分 割 的 精度 以 提高 组 织 对 比 度 。 第 二 ， 本 文 利 用 matlab 平台 实现 图 像 阴影 校正 算法 ， 按 表 1 的 参数 
校正 一 层 CBCT 图 像 需 2.1 秒 ， 其 中 Lo 范 数 平 滑 步 又 占用 1.5 Tb. Lo 范 数 的 不 可 解 性 导致 Lo 范 数 平 


滑 步骤 耗 时 ， 半 二 次 分 裂 的 特殊 交 蔡 优化 策略 增加 了 算法 的 时 间 复杂 度 ， 算 法 速度 仍 有 待 提高 。 


bem 


本 研究 提出 了 Lo 范 数 平滑 与 图 像 分 割 相 结合 的 CT 图 像 阴影 校正 方法 ， 该 方法 有 效 地 改善 了 CBCT 
的 图 像 质 量 。 在 均匀 模 体 实验 中 , 校正 后 图 像 ROI 的 CT 值 误 差 由 77 HU E 2 HU。 整体 图 像 的 SNU 
由 996 降 至 0.4%。 在 病人 头 部 试验 中 ， 校 正 后 图 像 ROI 的 CT 值 误 差 由 200 HU 降 至 13HU。 整 体 图 
像 的 非 SNU 由 1496 EZ 1%。 与 现 有 的 阴影 校正 方法 不 一 样 ， 本 方法 无 需 先 验 信息 ， 完 全 兼容 现代 放 
疗 加 速 器 的 图 像 引导 系统 ， 不 用 改变 其 他 硬件 和 扫描 协议 ， 不 增加 扫描 时 间 和 扫描 剂量 ， 具 有 一 定 的 临床 
应 用 价值 。 
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